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Abstract. - The scaling of the thermoremanent magnetization and of the dissipative part of 
the non-equilibrium magnetic susceptibility is analysed as a function of the waiting-time s for a 
simple ferromagnet undergoing phase-ordering kinetics after a quench into the ferromagnetically 
ordered phase. Their scaling forms describe the cross-over between two power-law regimes 
governed by the non-equilibrium exponents a and Xr/z, respectively. A relation between a, 
the dynamical exponent z and the equilibrium exponent n is derived from scaling arguments. 
Explicit tests in the Glauber-Ising model and the kinetic spherical model are presented. 
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The study of non-equilibrium critical phenomena originated from studies of the dynamic 
behaviour of glassy systems, but has also produced new insights into the behaviour of the 
conceptually simpler ferromagnetic systems, see for reviews. Since the latter are 

considerably more tractable, it might be hoped that insights gained from studying the kinetics 
of simple ferromagnets could provide clues for the comprehension of the former. In this letter, 
we consider simple ferromagnetic spin systems, quenched from a disordered initial state to 
below its critical temperature T c > 0. We shall work with a non-conserved order parameter 
throughout. 

Then the system undergoes phase ordering, that is domains of a time-dependent typical 
size L(t) ~ t x l z form and grow, where z is the dynamical exponent [2]- The resulting slow 
evolution of the system is more fully revealed through the study of two-time quantities, such 
as the two-time autocorrelation function C(t, s) and the autoresponse function R(t, s) 



C(t,s) = (<Kt)<t>(s)) , R(t,s)= 



(1) 

h=0 
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where <j) is the order parameter, h the conjugate magnetic field, t is called the observation 
time and s the waiting time. Ageing occurs in the regime when s and r = t — s > are 
simultaneously much larger than any microscopic time scale r m i cro . In many systems, one 
finds in the ageing regime a scaling behaviour, see 0E] 

C(t, s) = s- b f c (t/s) , R(t, s) = s-^fnit/s) (2) 

where a and b are non-equilibrium exponents and the scaling functions behave for large argu- 
ments x = t/s ^> 1 asymptotically as 

fc(x) ~ x~^' , f R (x) ~ z-W* (3) 

where Ac and A/j are the autocorrelation ISJE1 an( i autoresponse [7] exponents, respectively. 
For the usually studied case of an initial state without long-range correlations, it is gener- 
ally accepted that Ac = Ar = A. Combinations of rigorous results and of heuristic scaling 
arguments were used to derive the bounds d/2 < A < d [SlIHj- 

It has been proposed recently that one might be able to go beyond mere dynamical scale 
invariance as expressed by eq. J2J to a group of local scale transformations related to conformal 
transformations in time [§]. If that hypothesis applies, the form of the scaling function 

f R (x) = r x 1 + a - x */ z (x-l)- 1 - a (4) 

is completely fixed (r is a normalization constant). Eq. has been confirmed in several 
models, especially, through extensive simulations, in the 2D and 3D Glauber-Ising models |1U| . 

We are interested in the value of the non-equilibrium exponent a. In quite a few models, 
values of a were obtained, see [3] for a review, but no clear picture has yet emerged. In the 
exactly solved ID Glauber-Ising model, one has a = |1 11112) . However, the value of a in the 
2D and 3D kinetic Glauber-Ising models and the kinetic spherical models has been debated 
recently. In the 2D and 3D Glauber-Ising model, analytical ^Sj and numerical ^Oj results 
indicate a = 1/2. In the exactly solvable spherical model, one reads off a = d/2 — 1 for all 
spatial dimensions d > 2 from the exact result for R(t, s) )14ll5 16, 7, I^J. 

An alternative route has been followed by Corberi, Lippiello and Zannetti [17118119] . They 
studied the two-time correlation function C(t, s) and the susceptibility function s) in the 
O(n) vector model. Separating the two-time autocorrelation function C(t,s) = C s t(t — s) + 
Cage(i/ S ) m t° a stationary part and an ageing part and similarly for the ZFC susceptibility 
x(t,s) — Xst(t — s) + Xage^s), where Xst is defined such that it satisfies the fluctuation- 
dissipation theorem together with C cq , they extract an exponent a from the scaling behaviour 
Xage(^, s) = s~ a x(t/s) and claim that a = a. This was backed up by the assertion that the 
scaling of Xagc (t,s) were anomalous [E]. The same procedure was applied to simulational 
data of the ZFC susceptibility from the 2D and 31? Glauber-Ising model. If a n is the value 
of d in the 0(n)-model, they propose |17U18ll*H?| 

_ f (d-l)/4 ;d<3 / (d-2)/2 ;2<d<4 

ai - { 1/2 ; d > 3 ' 000 " \ 1 ; d > 4 (5j 

for the Glauber-Ising (n = 1) and the spherical model (n = oo), respectively. However, since 
floe 7^ d/2 — 1 = a for the spherical model, although even Corberi et al. do recover a = d/2 — 1 
from R(t,s), see eq. (83) in |18|. their claim of an anomalous scaling of x(t> s ) as raised 
in [TTnHl?j appears to be ill-founded. 

Our results in this letter are as follows: one should distinguish between those systems 
(called here class 5) with a short-ranged (exponential) decay of the equilibrium connected 
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spin-spin correlator in the ordered phase and those (called here class L) for which there 
remain long-range (algebraic) correlations in the ordered phase. The Glauber-Ising model in 
d > 1 belongs to class S while the spherical model and the 2D XY-model are members of class 
L |2I]. Then 

1/z ; class S , . 

{d-2 + rj)/z ; class L ^ ' 

where 77 is a well-known equilibrium critical exponent. Since for the Glauber-Ising model 
z = 2 2 j, one recovers the well-known a = 1/2 for class S in all dimensions d > 1. These 
results will be derived by comparing the scaling form of the dissipative part x" °f the non- 
equilibrium susceptibility (see eq. below) with heuristic expectations which for class S 
amount to the generally accepted idea that the ageing comes from the movement of the 
domain walls which separate the ordered domains m^lEI- I n addition, we shall derive a 
more complete scaling form for the thermoremanent magnetization Mtrm(^ s) m the limit of 
large waiting times s 1 

M TKM (t,s)/h= f duR{t,u) = s- a f M (t/s)+s- XR/z g M (t/s) (7) 



(i 



provided the system was initially prepared at infinite temperature. Eq. Q implies that the 
behaviour of Mtrm will be one of a cross-over between two distinct regimes I n 
practice, the cross-over region may well be large. If local scale invariance 9 holds, the scaling 
functions can be found from eq. iQJ and read explicitly (r^i are non-universal constants) 



f M {x)=r x- XR / z 2 F 1 (l + a,^-a;^-a + l;x- 1 ^ , g M {x) « n aT A «/« 



(8) 



The new and explicit forms of these scaling functions allow for a considerably more precise 
test of the cross-over scaling of Mtrm than had been possible before. In particular, having 
fixed the parameters r ,i for a given value of x, a precise prediction for the scaling of Mtrm 
for any other value of x, without any free parameter, is obtained. We shall use this idea to 
perform a new type of precision test on the scaling of Mtrm(^ s) for class S systems and in 
particular test for the value of the exponent a. 

We shall first derive the scaling forms of x" an d Mtrm as a function of the waiting time 
s and prove © by relating heuristic ideas on coarsening to scaling arguments. Then we shall 
describe numerical tests of (|7I8H in the spherical and Glauber-Ising models. 

Consider the time-dependent response of a ferromagnetic system to a sinusoidal magnetic 
field with frequency u. The imaginary part of the susceptibility is given by 00] 

X"(u;s)= [ du R(s,u) sin (uj(s - u)) (9) 
Jo 

which can be analysed using properties of the autoresponse function R(s,u). Here, the time 
difference r = s — u plays a central role and depending on its value an equilibrium or else an 
ageing behaviour is obtained. Specifically, it can be shown JS] that there is a time scale t p ~ 
such that < C < 1 (explicitly, £ = 4/ (d+2) for the spherical model [T2J on which the breaking 
of the fluctuation-dissipation relation occurs and furthermore i?(t, s) ~ R C q(t — s) for t — s<t p 
and R(t, s) = R sca i(t, s ) is given by © for i p < t — s. For u ~ s one measures the response 
with respect to a change in the initial conditions and R « Ri n i(t) ~ £~-W z [2]. Besides t p , 
and following E3E], we therefore introduce a third time scale t e such that s — t e = 0(1). 
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Changing variables and splitting the integral into three terms, one has 

x"(lo;s) = / dr R(s, s — r) sin lot 
Jo 

dr R(s, s — t) sin lot + / dr R(s, s — r) sin lot + / dr R{s, s — r) sin lot (10) 

ftp /-V s / 1 \ sinfwsw) /" s 
~ / dr R cq ( T ) sin lot + s~ a / du/ R V ^ + s~ Afl / z / drc sinter 

•/ Jt p /s \ L ~ V / \ i -— V ) Jt E 

where Co is a constant. In the third line, we used the asymptotic form of R(s, u) in the three 
regimes. This is justified if for sufficiently large values of s, the cross-over between these 
regimes is rapid. In a given system, this can be checked by looking for a rapid cross-over in 
the fluctuation-dissipation ratio (see e.g. 0120) between the initial equilibrium and the later 
ageing regime, see (TBIHS) . Now, in the limit s — > oo one has (i) t p — > oo, (ii) t p /s — » and 
(hi) t £ /s — > 1. We then find 

X >; s) = xi M + S - a X 2 (^s) + O ( S ~ A -/ 2 ) (11) 

For systems of class S, it is generally expected that the dynamics proceeds via domain 
growth, due to interface motion between well-ordered domains ^EIEIEI- Therefore, the 
susceptibility \ should decompose into a stationary part Xeq to which the spins in the bulk of 
the domains contribute and an ageing part Xag which comes only from the domain walls, with 
a length scale given by the typical distance L(s). The imaginary part of the susceptibility 
should read [D[Tl|[TSj 

x"{lo- S )= X ^{lo) + L{s)-\: s {lo S ) (12) 

Comparing with the general scaling form (|llfl . we see that the first term corresponds to the 
equilibrium term in l|12|) . the second to the ageing contribution and the third is a correction 
term. From the scaling L(s) ~ s 1 /* of the domains one indeed recovers the first part of 

For systems of class L the long-range fluctuations turn the ordered domains into fractal 
objects and we can no longer assume the existence of simple and well-separated domain walls. 
Rather, we anticipate instead of l|12f> the form 

X >; s) = x^H + L{ s )-^n) xlg{ujs) (13) 

(up to possible logarithmic factors, which arise e.g. in the 2D XY model [2])- The ageing 
part should come from the moving fractal domain boundaries. The simplest way to describe 
the time-dependence of this is to assume that all length scales are measured in terms of the 
domain size L(t). If that is so, the spin-spin correlator should read C(r,t) « C qs (r, L{t)) ^ 
r -( d - 2 +'n) c ( r J ' L{i)) in terms of a quasi-stationary correlator C qs and where c is some scaling 
function. The time-dependent contribution to the susceptibility per volume V should then 
come from Xqs ~ V ~ 1 J v dr C qs (r, L(t)) ~ L^)^^ 2 ^ and this leads to the desired result 
()13J) . The rest of the analysis then goes through as before and the second part of JHJ follows. 
Of course, we must assume that the scaling functions X& g ( x ) have no singularity as x — > oo. 

Eq. JBJ) is confirmed for class L by all known results for simple ferromagnets. First, for 
the short-ranged spherical model either without [14II15IIT51 I5] or else with |7] initial long- 
range correlations one has z — 2, 77 = 0. From the exact results for R(t, s) one recovers indeed 
a = d/2 — 1 as expected from JjjjJ. Second, for the spherical model with long-range interactions 
of the form J(r) ~ r~ d ~ a and < a < min(ci, 2) one has z = a, 77 = 2 — a and a = d/a — l |2*2^ . 
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Fig. 1 - Scaling of the thermoremanent magnetization Mtrm(£, s) (full curve) at T = 2 as a function 
of the waiting time s in the spherical model for x = t/s = 5. (a) d — 3 and So = 0, (b) d = 3 and 
So = 0.5, (c) d = 3.5 and So = 0.5. The dashed and dotted lines show the leading contributions 
~ s~ a and ~ s~ Xr ^ z , respectively, whereas the dash-dotted lines give their sum. 

as it should be. In these models, 6 = (see eq. |J2)). Third, in the 2D XY-model one has 
z = 2 and a = b = rj(T)/2 for both fully ordered and fully disordered initial conditions as 
shown analytically and through simulations |21j . where rj(T) depends on temperature. Fourth, 
the defining condition for class L is also satisfied for non-equilibrium critical dynamics. There, 
it is known that a = b = 2j3/vz = (d—2 + rf)/z, see [5], in agreement with JJJJ. Finally, for the 
ID Glauber-Ising model ageing occurs only at the critical point T c = 0. The known exponent 
r\ = 1 yields a = 0, as it should be 

In the same way, we now analyse the scaling of the thermoremanent magnetization. We 
stress that no use is made of eq. at this point. Introducing the same splitting of the 
integral, we find from the definition J7J) the scaling form Mtrm(£, s)/h = po + s~ a fAi(t/s) + 
s~ x ii/ z g M (t/s) where we used the asymptotic forms of R of the three regimes. Then Jm(x) 
is related to Jr(x) |5j and the last term (and consequently guix), see eq. (jSJl) was roughly 
estimated from the mean-value theorem. The term po is related to the initial state of the 
system. In particular, po = for an initially fully disordered state. While this expression for 
Mtbm depends only on the scaling J5J), from local scale invariance as introduced and tested 
in [HlllOj we have eq. (@J and recover JHJ. 

We are now ready for numerical tests of Q and JSJ. We begin with the kinetic spherical 
model, formulated in terms of real spin variables S x (t) subject to the (mean) spherical con- 
straint J2 X (Sx(t) 2 ) = A/" where J\f is the number of sites (see [22] for a careful discussion on 
this point). The kinetics is given through a Langevin equation 

= E S v^ (2d + m)S*(t)+V*(t) (W) 

where the sum over y extends over the nearest neighbours of x. The thermal noise is assumed 
to satisfy {rj x {t)-q y {t')) = 2TS x y S(t — t'). From the spherical constraint the Lagrange mul- 
tiplier i(t) can be found through the solution of a Volterra integral equation, and two-time 
correlators and response functions are readily found, see [mrH)ll23| for details. 
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Fig. 2 - Scaling of the thermoremanent magnetization Mtrm(^, s) as a function of s for several values 
of x — t/s (squares: x = 5, circles: x — 7, triangles: x — 9), as obtained in the 2D Glauber-Ising 
model at T — 1.5. The black curves are calculated from equations J7J and JHJ with Xr — 1.26, 
see |5Hrn) : (a) a = 1/2 and (b) a = 1/4. 



In figure ^ we show the thermoremanent magnetization Mtrm obtained from numerically 
integrating the Langevin equation, at a temperature T < T c (recall that for d — 3, T c ~ 3.96 
and for d = 3.5, T c ~ 5.27). As an initial state, we used uncorrelated spins with either a mean 
magnetization So — (figure QJi) or else So = 0.5 (figure^ andQJ;). Then the exponents 
z = 2, a = d/2 — 1 and X R = d/2 for S — or X R = d for S = 0.5, respectively, are 
expected \7\. The leading term ~ s~ a and the sum of the two leading contributions to Mtrm 
according to eq. (J7J) are shown and are compared with the numerical solution of the Langevin 
equation - For So = 0.5, the cross-over between the two regimes is evident. On the other 
hand, for So = 0, we observe that although the slope of In Mtrm versus Ins appears to be 
fairly constant, the second term expected from eq. (0 produces a sizable correction. We shall 
find a similar behaviour for the 2D Glauber-Ising model below. Finally, we see that down to 
values of s as small as s ~ 1, the time-dependence of Mtrm is well described by eq. J7J). We 
observed a similar behaviour, and in quantitative agreement with the scaling functions JHJl, 
for other values of x. 

We now come to our main objective: deciding between a = 1/2 and a = 1/4 through a 
test of (0) in the 2D Glauber-Ising model. This model describes the evolution of the spin 
variables Sj = ±1 on the sites j of a hypercubic lattice, realized through heat-bath Glauber 
dynamics. We use a temperature T = 1.5 < T c and a spatially random field h — ±0.05 and 
measure Mtrm(*) s ) m the standard way I2H1- Starting from a totally disordered state, we 
measure the evolution of Mtrm(£, s) as a function of s with x = t/s fixed. The non- universal 
parameters r^i were fixed using the data for x = 7 for either a = 1/2 (tq = 1.76, n = —1.84) 
or a = 1/4 (ro = 0.22, r% = 0.09), respectively. Since ro,i are independent of x, we therefore 
obtain a prediction on how Mtrm(^s, s) should scale for any other value of x. This is shown 
in figure|3 While for a — 1/2, the Monte Carlo data nicely agree with the cross-over scaling as 
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predicted from eqs. I|7I8|) . it is also clear from figure 03 that the hypothesis a = 1/4 |18lll9j is 
incompatible with the data. The same conclusion, namely the invalidity of the value a = 1/4 
in the 2D Glauber-Ising model, has also been reached using simulational data at T = with 
waiting times up to s — 5600 27 . 

Summarizing, we have reanalysed the scaling of the magnetic linear response in simple fer- 
romagnets undergoing phase-ordering kinetics. The cross-over scaling of the thermoremanent 
magnetiztion eq. J7J) (and also the dynamic susceptibility) is predicted exactly for any dynamic 
universality class from local scale invariance in terms of the universal exponents Xr/z and a. 
We have confirmed eqs. © and (|7I8|) in the Glauber-Ising and spherical models. This provides 
evidence against recent suggestions of an anomalous scaling of the dynamic suceptibilities in 
phase-ordering kinetics. 
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